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Abstract 

Background: In this study, we developed a pharmacokinetic (PK)- pharmacodynamic (PD) model of a new sustained 
release formulation of interferon-a-2a (SR-IFN-a) using the blood concentration of IFN-a and neopterin in order to 
quantify the magnitude and saturation of neopterin production over time in healthy volunteers. The SR-IFN-a in this 
study is a solid microparticular formulation manufactured by spray drying of a feeding solution containing IFN-a, a 
biocompatible polymer (polyethylene glycol) and sodium hyaluronate. 

Methods: The full PK and PD (neopterin concentration) datasets from 24 healthy subjects obtained after single doses 
of 9, 18, 27 and 36 MIU of subcutaneous SR-IFN-a were used to build the mixed-effect model using NONMEM (version 
7.2) with the G FORTRAN compiler. 

Results: A one-compartment model with first-order elimination and a mixture of zero- and first-order absorption was 
chosen to describe the PK of SR-IFN-a. The time-concentration profile of neopterin, the PD marker, was described by a 
turnover model combined with a single transit compartment. The saturable pattern of the neopterin response blurring 
the dose-response relationship of SR-IFN-a was addressed by introducing the concept of the EC 50 increasing over time. 

Conclusions: The PK-PD model of SR-IFN-a developed in this study has presented a quantitative tool to assess the 
time-course of a saturable neopterin response in humans. 
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Background 

Interferons are produced by peripheral blood leukocytes, 
fibroblasts, and activated T and NK cells in response to 
viral infection or other inducers including double- 
stranded RNA, lipopolysaccharide, micro-organisms, or 
endotoxins [1-3]. Based on immunological and physico- 
chemical differences, human interferons are divided into 
a-, (3-, and y-interferon families, with numerous subtypes 
within each interferon family [4]. Because of its antiviral, 
antiproliferative and immunomodulating properties, re- 
combinant interferon-a (IFN-a) has been used as a 
treatment for various diseases [5]. However, in the case 
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of chronic hepatitis C, monotherapy with IFN-a has 
been persistently effective in only a small percentage of 
patients [6]. This low response rate is thought to be due 
to HCV genotype variation and/or the quite short half- 
life of IFN-a [7,8]. Combination therapy with other anti- 
viral agents, such as ribavirin, is therefore recommended 
[9]. Frequent administration (3 times weekly) of IFN-a 
has been considered to be an additional cause for thera- 
peutic failure of interferon due to the fact that frequent 
administration accelerates the formation of neutralizing 
antibodies and causes other adverse effects due to large 
variations in peak-to-trough plasma drug concentrations 
[10,11]. Hence, a long-acting formulation has been de- 
veloped for IFN-a which can extend its effects to weeks 
and months. 

Unlike small molecule drugs, the poor stability of pro- 
tein drugs has been a hurdle to the development of long 
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acting formulations [12]. Technologies for producing 
long-acting formulations of protein drugs can be roughly 
classified into chemical modification such as pegylation 
[13] or formulation changes allowing delayed release 
from depot sites. The new sustained release formulation 
of IFN-a-2a (SR-IFN-a, LG Life Sciences) used in our 
report is a solid microparticular formulation manufactured 
by the spray drying of a feeding solution containing IFN-a, 
a biocompatible polymer (polyethylene glycol) and sodium 
hyaluronate [14]. 

This report is based upon a first-in-human, single as- 
cending dose trial of SR-IFN-a in healthy volunteers 
where the IFN-a concentration profile showed an 
extended release pattern through the doses studied. 
Blood neopterin concentrations were also measured as a 
biomarker demonstrating the activity of IFN-a in this 
first-in-human study. Neopterin is a soluble immune ac- 
tivation marker released from monocytes and macro- 
phages by IFN-a [15]. Looking into the relationship 
between the exposure to IFN-a and neopterin in the 
healthy subjects' data, we found that the magnitude of 
neopterin concentration changes was not clearly corre- 
lated with the dose of IFN-a. Although similar phenom- 
ena have been observed in animal experiments [16,17], 
this has never been reported in humans despite frequent 
clinical trials using neopterin as a biomarker. Because 
clear understanding of the relationship between exposure 
and response is one of the fundamental goals of early- 
phase exploratory clinical trials, we investigated the 
concentration-response relationship of IFN-a and neopterin 
in humans, which has never been elucidated. As results, 
we present a pharmacokinetic-pharmacodynamic (PK-PD) 
model that quantifies the peculiar time-course of 
neopterin responses to IFN-a in humans. 

Methods 

Inclusion and exclusion criteria 

Healthy volunteers aged 18 to 45 years with BMI ran- 
ging from 19 to 29 kg/m 2 , with no clinically relevant 
conditions identified based on medical history, physical 
examination, laboratory tests or electrocardiography 
(ECG), were eligible for inclusion. Subjects with any history 
that indicated a possible alteration in IFN-a metabolism 



or with hypersensitivity to IFN-a were excluded. The final 
study enrollment was 32 subjects (Table 1). 

Study design 

A randomized, double-blind, active controlled, dose es- 
calation phase I clinical study was conducted on 32 
healthy subjects in the clinical pharmacology unit of the 
Kendle International BV, located in Utrecht, Netherlands. 
Subjects were randomly allocated into four groups (eight 
subjects per group). Within each group, six were given 
SR-IFN-a (test formulation) and the other two were given 
3 MIU Roferon-A® (Roche, active comparator) via sub- 
cutaneous injection. The doses of SR-IFN-a allocated to 
groups 1, 2, 3 and 4 were 9 MIU, 18 MIU, 27 MIU, and 36 
MIU, respectively. 

The study was performed in compliance with the 
European Community rules of Good Clinical Practice 
(GCP), the International Conference on Harmonization 
(ICH) Tripartite Guidelines: Guideline for GCP, the current 
revision of the 'Declaration of Helsinki' (Edinburgh, 
amendment October 2000). The Stichting Therapeutische 
Evaluatie Geneesmiddelen (STEG), an independent ethics 
committee, approved the protocol before execution of the 
trial, and all participants gave written informed consent. 

Blood sampling 

For the population PK analysis, peripheral blood samples 
(5 mL each) were taken just prior to the injection, and 
0.75, 1.5, 3, 6, 8, 10, 12, 18, 24, 30, 36, 48, 60, 72, 96, 
120, 144, 168, and 192 hours after the injection. For the 
PD marker, the sampling times differed slightly: just 
prior to the injection and 3, 8, 12, 18, 24, 36, 48, 72, 96, 
120, 144, 168, 192, and 264 hours after the injection. 
The samples were collected in light-protective tubes and 
stored at < -70°C. 

Assay of plasma concentrations of IFN-a and neopterin 

IFN-a concentrations in the serum samples were ana- 
lyzed using a commercial Human IFN-a Multi-Subtype 
ELISA Kit (product # 41105) with a detection limit of 
12.5 pg/mL manufactured by Pestka Biomedical Laborator- 
ies, Inc. (Piscataway, NJ, USA). Neopterin concentrations 
in the serum samples were analyzed using a commercially- 



Table 1 Subject demographies' 1 



Group 


Control 


1 




2 


3 


4 


Dose 


3 MIU 


9 MIU 




18 MIU 


27 MIU 


36 MIU 


Number of Subjects 


8 


6 




6 


6 


6 


Sex 


Male 


Male 




Male 


Male 


Male 


Age (years) 


24.0(18 - 43) 


21.5 (19 ~ 


37) 


22(19 - 34) 


22 (18 - 42) 


24 (20 ~ 43) 


Height (cm) 


181.65 (165.5 ~ 189.0) 


179.05 (171.0 - 


- 191.5) 


183.5 (172.0 ~ 187.0) 


178.15 (168.0 ~ 191.5) 


181.75 (178.0 - 193.0) 


Weight (kg) 


73.85 (69.0 ~ 82.2) 


78.65 (60.3 - 


■ 92.4) 


76.15 (63.0 ~ 86.5) 


74.40 (53.8 - 92.6) 


74.40 (70.7 - 98.8) 



Continuous variables are shown as median (range). 
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available Neopterin ELISA method (REF 40-371-25012, 
GenWay Biotech, Inc., San Diego, CA, USA) with a detec- 
tion limit of 0.7 nmol/L and the specificity of about 
99.95%. 

Population PK-PD model 

Because the aim of this study was to develop a PK-PD 
model for SR-IFN-a, data from the active control group 
participants, who were given the immediate release IFN- 
a formulation (8 subjects), were not included in the ana- 
lysis (individual plots for PK-PD models are shown in 
Additional file 1). 

Mean plasma concentrations of IFN-a are shown in 
Figure 1, and non-compartmental analysis results of the 
PK of SR-IFN-a are summarized in Table 2. The popula- 
tion PK-PD analysis was performed using NONMEM 
(version 7.2, Icon Development Solutions, Ellicott City, 
MD, USA) with the GFORTRAN compiler. 

To find the model that best described the absorption 
profile, which showed double peaks in many subjects, 
first- and zero-order absorption models and their com- 
bined form, with or without lag time, were tested. Based 
on first-order elimination, one- and two-compartment 
distribution models were tested for the three absorption 
(first-order, zero-order and combined) processes. The 
Michaelis-Menten absorption and elimination models 
were also tested considering the potential saturable ab- 
sorption or elimination using the ADVAN subroutines 
(Table 3). 

The coupling of the PK model to the PD model was 
done in a sequential manner. The population PD model- 
ing was performed using the individual PK parameters 
estimated from the final PK model, which were added to 
the PD dataset. 



Turnover models, with or without transit compart- 
ments, were compared to find the most appropriate 
model that explained the delayed effect of IFN-a on 
neopterin concentrations. The turnover model was ini- 
tially selected over the effect compartment model based 
upon the well-known action of IFN-a stimulating the re- 
lease of neopterin, and the transit compartments were 
tested because their usefulness was reported in a previ- 
ous preclinical study [16]. 

The differential equations for the drug effect model 
(model structures are shown in Figure 2) were: 



dA(l) 
dt 



dA{2) 
dt 



dA(3) 
dt 



dA{4) 
dt 



-K a -A(l) 



■ K a -A(l)-K e -A(2) 



■ K in -E(C)-K tr 'A{3) 



■ K tr -A(3)-K out -A(4) 



where K a and K e are the absorption and elimination rate 
constants for IFN-a, respectively. K in is the production 
rate of neopterin, a zero-order constant, K tr is the first- 
order transition rate leaving the transit compartment, 
K out is the first-order rate constant for the elimination of 
neopterin, and E(C) is the effect as a function of the in- 
dividual predicted drug concentration, C. A(4), which is 
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Table 2 Non-compartmental analysis of interferon -a: Mean ± S.D. (CV%) 



Group 


Dose 


Formulation 


C m ax (pg/mL) 


Median T max (h) 


AUC, ast (ng.h/mL) 


Control 0 (N=8) 


3 MIU 


Roferon-A® 


43.27 ± 9.804 (22.66%) 


6 


0.9664 ± 0.9322 (96.46%) 


1 (N=6) 


9 MIU 


SR-IFN-a 


28.33 ± 9.656 (34.08%) 


18 


2.072 ±1.134 (54.73%) 


2 (N=6) 


18 MIU 


SR-IFN-a 


62.12 ± 15.93 (25.64%) 


24 


5.373 ± 1.382 (25.72%) 


3 (N=6) 


27 MIU 


SR-IFN-a 


65.73 ± 6.702 (10.20%) 


24 


5.544 ± 0.5509 (9.94%) 


4 (N=6) 


36 MIU 


SR-IFN-a 


80.31 ± 9.859 (12.28%) 


24 


7.151 ± 1.132 (15.83%) 



a There was not a separate control group, but data from a total of eight subjects (two subjects in each of the four groups) who received 3 MIU of Roferon-A { 
were shown. 



the amount existing in the 4th compartment, indicates 
the concentration of neopterin. 

The stimulatory function for the drug effect, E(C), was 
a sigmoid function: 



E{C) = 1 + 



EC 



50 



where E max is the maximum effect and EC 50 is the IFN- 
a concentration that produces 50% of the maximum 
effect. 

In our study, differences in the mean plasma neopterin 
concentrations between the groups receiving different 
doses of SR-IFN-a were not clearly discernible (Figure 3). 
To account for this phenomenon in our PD model, we 
incorporated the concept of time-dependent attenuation 
of the effect parameters [18], especially the increasing 
EC 50 over time, as shown in the following equation: 



EC 50 = ECB x (l + CA x (l-e 



-CBxTIME 



)) 



where ECB is the baseline of EC 50 and CA and CB are 
coefficients to describe the EC 50 increase to a certain 
level (ECB x CA) over time. This curve is the cumula- 
tive distribution function (CDF) of an exponential distri- 
bution to EC 5 o It has advantages in explaining the 
concave-shaped curves in relation to time. 

The first-order conditional estimation (FOCE) method 
with interaction was used throughout the model build- 
ing. Models were selected based upon a decrease in the 



objective function value (OFV) of more than 3.84 (P-value 
0.05 in an approximate j2 distribution) and improvement 
in the individual plots, as well as other scatterplots. 

A log normal distribution was assumed for inter- 
individual variability (n), and PK or PD parameters of 
the jth subject (Pj) were described as: 



Pj = TVP x exp(/7 ; .) 



where TVP represents the typical population value of 
PK-PD parameters, such as clearance (CL), volume of 
distribution (V), absorption rate constant (K a ), lag time 
(ALAG), the first-order elimination rate of serum 
neopterin (K ou t)> an d maximum stimulation effect (E max ). 
The inter-individual variability eta (n) for each PK-PD par- 
ameter was assumed to follow a Gaussian distribution 
with a mean 0 and a variance co 2 . Possible correlations be- 
tween the inter-individual variability were also evaluated. 

As for the residual error, the additive, proportional 
and combined forms were tested. An example of the 
combined error form is shown as follows: 



Yij = IPREDtj + 



where IPREDij is the individual predicted concentration, 
is the measured concentration of the j th individual at 
the i th sampling time, and is residual error. Residual 
errors (e) include intra-individual variability, assay 
error and model misspecification. They were also 



Table 3 PK model development process 

Step 



Model tested 0 



Objective function value 



1 One-compartment model with first-order absorption 2138.416 
One-compartment model with first-order absorption (with lag-time) 2130.609 

2 Two-compartment model with first-order absorption 2037.527 
Two-compartment model with first-order absorption (with lag-time) 1993.071 

Two-compartment model with saturable absorption (Michaelis Menten absorption) 2140.583 

3 Two-compartment model with zero-order absorption 1885.795 

4 One-compartment model with a mixture of zero and first-order absorption 1879.067 
One-compartment model with Michaelis Menten elimination and a mixture of zero- and first-order absorption 2163.942 



a The elimination process was assumed to follow first-order kinetics. 
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PK domain 



Depot (1) 



zero- \ 
order 1 ) V 



K 2) 



Central (2) 



PD domain 





k in 


Transit (3) 




K 


Neopterin(4) 



1 >0<time<D2 : zero-order input (F z -F-Dose) 

2 > From time ALAG : first-order input ((1-F z )-F-Dose) 

©: Stimulatory effect attenuated by time 

Figure 2 PK-PD model structure. The mixed absorption with one- 
compartment first-order elimination model (PK) and the turnover 
model with a single transit compartment (PD) are shown. 



assumed to follow a Gaussian distribution with a mean 0 
and variance a 2 . 

Covariate selection 

Age, height, weight and creatinine clearance were 
screened as potential covariates of the parameters using 
Generalized Additive Modeling (GAM) implemented by 
Xpose version 4.2.3. 

In the forward selection of covariates, variables that 
decreased the OFV by more than 3.84 (P < 0.05) and im- 
proved the inter-individual variability (omega value de- 
crease) were selected. Covariates that did not increase 
the OFV more than 3.84 (P < 0.05) in the backward 
elimination step were removed from the model. 



Model evaluation 

The 95% confidence intervals (CIs) for mean population 
PK and PD parameters were determined using a re- 
sampling technique based on the bootstrap method. 
One thousand re-sampled datasets were collected and 
their parameters were estimated using our population 
models. The models were also evaluated by visual pre- 
dictive checks (VPC) using 1,000 simulated datasets. 

Results 

Population PK-PD model 

The results of non-compartmental PK analysis showing 
the trend of PK linearity of SR-IFN-a are briefly summa- 
rized in Table 2. A one-compartment model with first- 
order elimination and a mixture of zero- and first-order 
absorption best described the PK of SR-IFN-a. The de- 
layed pattern of the time-concentration profile of 
neopterin was well described by the turnover model with 
a single transit compartment. There was no significant 
covariate. The structure of the final PK-PD model is 
shown in Figure 2. The population PK-PD parameter es- 
timates, with corresponding standard error (SE) values, 
are summarized in Table 4. Basic goodness-of-fit plots 
are presented in Figure 4. Predicted concentration-time 
profiles of IFN-a and observed data from representative 
individuals are shown in Figure 5. 

Model evaluation 

The 95% confidence intervals (CIs) for population PK- 
PD parameter estimates, determined using the bootstrap 
re-sampling method, are shown in Table 4. VPC plots 
simulated concentrations of 1,000 virtual datasets 



Neopterin 



Roferon A, 3 MIU 
SR-IFN-A, 9 MIU 
SR-IFN-A, 18 MIU 
SR-IFN-A, 27 MIU 
SR-IFN-A, 36 MIU 




Figure 3 Plasma neopterin concentration changes in response to interferon-a (Mean and S.D.). Plasma neopterin concentrations after 
injection of Roferon-A® or SR-IFN-a show little difference between dose groups. 
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Table 4 Final estimates of population PK-PD parameters 

Parameter Description Estimate %RSE Bootstrap median (95% Cl)° 

Pharmacokinetics 



CL/F (L/h) 


Apparent clearance 


12.2 


7.39 


12.5 (4.9-14.9) 


V/F (L) 


Apparent volume of distribution 


691 


6.54 


712 (324-980) 


D2 (h) 


Duration of zero-order absorption 


20.2 


7.08 


19.9 (16.6-23.0) 


KA (h~ 1 ) 


Absorption rate constant of first-order absorption 


0.00653 


16.23 


0.00721 (0.0010-0.1715) 


ALAG (h) 


Lag time to the initiation of first-order absorption 


85.7 


3.92 


88.1 (80.9-108.0) 


RF b 


Alternate variable for relative fraction absorbed in zero-order absorption process 


0.185 


55.68 


0.24 (-1.17-0.94) 


OOcL (%) 


Between subject variability of CL 


26.1 


35.1 


24.5 (14.4-32.7) 


oo v (%) 


Between subject variability of V 


23.8 


58.0 


23.7 (0.4-37.6) 


U) D2 (%) 


Between subject variability of D2 


35.7 


27.1 


34.8 (23.0-49.4) 


OOrf (%) 


Between subject variability of RF 


34.7 


64.1 


32.1 (0.3-54.2) 


ook a (%) 


Between subject variability of KA 


77.0 


37.1 


59.5 (0.7-95.4) 


Oadd (pg/mL) 


Additive error 


3.92 


12.65 


3.83 (-3.36-4.80) 


Oprop (%) 


Proportional error 


7.8 


26.06 


7.0 (-8.8-11.9) 


Pharmacodynamics 








BASE (nmol/L) 


Baseline of neopterin 


5.85 


4.17 


5.88 (5.39-6.34) 


KOUT (h" 1 ) 


First-order elimination rate of serum neopterin 


0.0311 


17.20 


0.03 (0.02-0.04) 


EMAX 


Maximum stimulatory effect 


16.1 


53.19 


1 7.45 (7.06-62.42) 


GA 


Hill coefficient (y) 


1.24 


11.85 


1.27 (1.01-2.52) 


CA C 


Coefficient 


405 


115.80 


616.5 (27-13115) 


CB C 


Coefficient 


0.0068 


27.73 


0.0064 (0.0014-0.0119) 


ECB C 


Baseline of EC 50 


2.17 


152.07 


1 .66 (0.05-74.67) 


M^ (h) 


Mean transit time 


14.6 


11.37 


14.5 (9.82-19.2) 


U) B ASE (%) 


Between subject variability of BASE 


13.85 


51.84 


12.63 (3.91-20.79) 


OOCB (%) 


Between subject variability of CB 


57.31 


43.31 


46.21 (0.32-72.61) 


03 G A (%) 


Between subject variability of GA 


13.51 


55.80 


13 (0.32-36.49) 


ooecb (%) 


Between subject variability of ECB 


21.31 


172.97 


21.39 (0.45-55.86) 


03 M ^ (%) 


Between subject variability of M^ 


13.36 


94.92 


10.75 (0.55-24.21) 


o add (nmol/L) 


Additive error 


1.14 


11.05 


1.11 (0.95-1.37) 



a 95% CIs obtained from estimation of 1000 bootstrap-resampled datasets; ^Fraction absorbed in zero-order absorption, F z = e RF /(1+ e RF ); Fraction absorbed in first- 
order absorption: 1-F Z ; C EC50 =ECB*(1+CA*(1-e~ CB * Time )); RSE Relative standard error. 



(nsub=1000 in the $SIMULATION block, 24,000 virtual 
patients) from the final model. The results from the 
VPC showed that the PK-PD model gave acceptable pre- 
dictive performance. Curves for the 12.5th, 50th and 
87.5th percentiles of concentrations were overlaid on the 
observed concentrations (Figure 6). 

Discussion 

In this study, we presented results of PK-PD modeling of 
the time-concentration profiles of IFN-a and neopterin 
after administration of SR-IFN-a in healthy subjects. 

For PK characteristics, the combined absorption 
model was successful in describing the double peak 
phenomenon. A zero-order absorption model appropri- 
ately described the initial increase in concentration, as 



measured by relatively frequent sampling, before 
reaching the maximum concentration. Subsequent sec- 
ond peaks (observed around 100 h after injection) were, 
however, best described by a first-order absorption 
model. One possible explanation for this double-peak 
phenomenon is in the method of SR-IFN-a administra- 
tion: SR-IFN-a should be mixed with medium chain tri- 
glycerides (MCT) right before subcutaneous injection. 
Micro-droplets with various sizes might be formed in 
this mixing step, and their absorption rates may differ by 
droplet diameters. Despite differences in formulation, 
the CL (12.2 L/h) and V (691 L) in our report were not 
much different from those reported in a previous study 
for IFN-a (7-8 L/h and 700-850 L, respectively) in 
healthy subjects [19]. 
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(A) IFN-a 





0 20 40 60 80 100 
Population Prediction (pg/mL) 



20 40 60 80 
Individual Prediction (pg/mL) 



20 40 60 80 100 
Individual Prediction (pg/mL) 




Time after Dose (h) 



(B) Neopterin 





Population Prediction (nmol/L) 



Individual Prediction (nmol/L) 





Individual Prediction (nmol/L) 



50 100 150 200 250 
Time after Dose (h) 



Figure 4 Basic goodness of fit plots for PK and PD models. (A) and (B) are goodness of fit plots for PK (IFN-a) and PD (neopterin) models, 
respectively. Black line, line of identity; gray line, LOESS (locally weighted regression) smooth line. The encircled dots in the PD plots represent 
one outlier (ID No. 2) whose IFN-a concentrations were very low, but whose neopterin concentrations were rather higher. 



As for the PD model, the addition of a transit com- 
partment before the neopterin compartment gave better 
outcome when compared with a simple turnover model 
The transit compartment for neopterin was used to 
model data from monkeys [16], and we found that it is 



also useful in a human PK-PD model in this study. As 
there are a few mechanistic models that tried to explain 
little PD differences between dose groups (saturation of 
responses) [16,20,21], we tested them for our neopterin 
data in the preliminary PD model development step; 
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Group2(18MIU) 



Group 3 (27 MIU) 



Group 4 (36 MIU) 



Medians of All dose levels 




150 200 



Time after Dose (h) 



50 100 150 200 

Time after Dose (h) 



50 100 150 200 

Time after Dose (h) 



50 100 150 200 

Time after Dose (h) 



50 100 150 200 

Time after Dose (h) 



° Observed value 

87.5% prediction interval 

Median value of prediction 

12.5% prediction interval 



0 50 100 150 200 250 0 50 100 150 200 250 0 50 100 150 200 250 0 50 100 150 200 250 0 50 100 150 200 250 

Time after Dose (h) Time after Dose (h) Time after Dose (h) Time after Dose (h) Time after Dose (h) 

Figure 6 Visual predictive checks of PK (Top) and PD (Bottom) models. One thousand datasets (24,000 virtual subjects) were simulated 

using the final PK and PD parameter estimates. The simulated median (solid lines) and 75% prediction intervals (broken lines) were overlaid with 

observed data. The two rightmost panels presenting medians of all dose groups show dose-linearity in PK, but not in PD. 
k J 



however, none of them were successfully converged by 
NONMEM. Neither the precursor turnover model [20] 
that explained tolerance with depletion of precursor 
molecules, nor the turnover feedback model [21] that 
explained tolerance with negative feedback via a moder- 
ator compartment provided acceptable parameter esti- 
mates, and the basic goodness of fit plots were even 
worse than those for the descriptive model of increasing 
EC 50 . The inhibitory feedback model in monkeys [16] 
also showed similar problems. Thus, we had to use the 
concept of time-dependent attenuation of EC 50 that did 
not include mechanistic reasoning to describe our 
neopterin response. 

There are conflicting reports on saturation of neopterin 
by IFN-a in patients. A report on patients with hairy cell 
leukemia showed that neopterin responses and clinical 
efficacies after low doses (0.5-0.8 MIU/ day) were similar 
to those after a conventional dose (3 MlU/day) in a 6- 
month clinical trial [22]. In another clinical trial for a con- 
trolled release formulation of IFN-a, increases in the mean 
AUC of neopterin were marginal among the doses tested 
(4.09 mM-h after 20 \ig IFN-a and 6.61 mM-h after 320 \ig 
IFN-a) [23]. The forms of IFN-a used in those studies 
[22,23] were non-pegylated forms, like the SR-IFN-a used 
in this report; however, in phase I clinical trials of 
pegylated IFN-a, neopterin responses were well correlated 
with the doses used [24,25]. Such a discrepancy in 
neopterin responses between non-pegylated and pegylated 
IFN-a suggests that the polyethylene glycol tail attached 
to IFN-a changed its PD parameters related to neopterin 
production. To the best of our knowledge, although 
pegylated formulations have long been used, it has never 



been reported that their PD profile may be different from 
that of non-pegylated forms. 

Because neopterin is a frequently-used marker of cell- 
mediated immunity, it can be used to monitor the de- 
gree of immune activation in various clinical conditions, 
including infections, autoimmune diseases, malignancies, 
and other conditions [26] . Neopterin is also known to me- 
diate the cytotoxic action of activated macrophages and 
dendritic cells via interactions with reactive oxygen species 
[27]. Thus, saturation of the neopterin response suggests 
that the magnitude of cytotoxic action mediated by 
neopterin may be similar regardless of the doses of IFN-a. 

Conclusions 

We developed a human PK-PD model revealing the sat- 
urable neopterin response to IFN-a for the first time. 
Our model suggests that the magnitude of cytotoxic ac- 
tion mediated by neopterin may be similar regardless of 
the doses of non-pegylated IFN-a. 
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